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Abstract : 

A  formula  is  derived  for  determining  dynamic  stress  intensity  factors 
directly  from  crack  mouth  opening  displacements  in  dynamic  tear  test  specimen. 
The  results  obtained  by  the  present  estimation  method  for  stationary  as  well 
as  propagating  cracks  agree  excellently  with  those  directly  obtained  through 
a  highly  accurate  moving-singularity  finite  element  method.  The  present 
method  can  also  be  applied  for  other  types  of  specimen  which  have  a  relatively 
short  edge  crack  without  any  loading  on  the  crack  surface.  The  present  simple 
estimation  method  should  be  of  great  value  in  the  experimental  measurement  of 
dynamic  stress- Intensity  factors  for  propagating  cracks  in  (opaque)  structural 
steel  dynamic  tear  test  specimens. 
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i' T »>ai >  .  >timation  Method  : 

There  is  a  major  difficulty  to  determine  experimentally  the  dynamic  stress 
intensity  factors  for  a  propagating  crack  in  a  non-transparent  (metallic)  frac¬ 
ture  specimen  under  the  conditions  of  impact  loading  as  well  as  quasi-static 
loading.  For  a  transparent  (photoelastic)  specimen,  dynamic  stress  intensity 
factors  can  be  measured  by  means  of  the  optical  method  such  as  the  photoelastic 
method  [1]  or  the  shadow  optical  method  of  caustics  [2].  On  the  other  hand, 
for  a  non-transparent  specimen,  so  far,  only  one  attempt  has  been  made  for 
determining  dynamic  stress  intensity  factors  directly  from  the  shadow  optical 
method  of  caustics  [3],  This  was  done  by  measuring  the  light  reflected  from 
the  mirrored  surface  of  the  high  strength  steel  specimen  with  a  dynamically 
propagating  crack  under  conditions  of  quasi-static  loading.  However,  the 
experimental  results  show  large  oscillations  of  stress  intensity  factor  with 
high  frequencies.  These  oscillations  are  attributed  to  the  high  frequency 
stress  waves  which  may  exist  only  on  the  surface  of  the  specimen.  This 
situation  may  become  more  critical  in  the  impact  test  specimen. 

In  the  present  note,  a  simple  method  is  developed  for  determining 
dynamic  stress  intensity  factors  directly  from  the  measurement  of  crack  opening 
displacement  at  the  notch  mouth  in  dynamic  tear  testing.  A  formula  for  the 
relation  between  dynamic  stress  intensity  factor  and  crack  opening  displace¬ 
ment  at  the  notch  mouth  is  derived  by  using  the  relation  in  the  static  case. 
This  possibility  was  found  from  a  critical  examination  of  the.  results  for 
several  numerical  fracture  simulations  mad o  in  Ref.  [41.  The  moving-s ingu- 
laritv  dement  procedure  used  in  Ref.  [4]  give.1;  a  direct  evaluation  of  the 
ivna.r.ic  'tress  intensity  factors  for  a  propagating  as  welL  as  stationary 
c  r  ac  iv . 

The  notations,  and  dimensions  of  the  dymunic  tear  rest  specimen,  wire  shown 
ia  'ig.  L.  fhe  specimen  geometry  and  material  properties  correspond  to  those 
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reported  in  Ref.  [5].  In  the  finite  element  analysis  [4],  the  following  initial 

▲ 

conditions  were  used:  at  time  t=0,  velocity  u^=6.88  m/sec.  The  tup  displace¬ 
ment  is  calculated  by  uT=uTt#  The  displacements  at  the  supports  are  fixed  for 
all  the  times.  In  Ref.  [A],  the  influence  of  the  loss  of  contact  of  the  speci¬ 
men  at  various  times  with  either  the  tup  or  the  supports  has  also  been  investi¬ 
gated.  However*  we  shall  focus  our  attention  to  the  case  when  the  specimen 
is  in  contact  with  the  tup  and  the  supports.  A  plane-strain  condition  is  invoked 
in  the  two-dimensional  analysis. 

The  crack  propagation  history  is  assumed  as  shown  in  Fig.  2.  This  crack 
propagation  history  was  used  as  input  data  for  the  "generation  phase"  fracture 
simulation.  As  shown  in  Fig.  2*  the  crack  is  stationary  (C=0)  until  t-95  usee, 
during  the  period  of  95  to  146  psec  the  crack  propagates  with  the  speed  of 
37.5  m/sec,  and  after  146  usee  the  crack  speed  becomes  95  m/sec. 

The  crack  opening  profiles  as  well  as  the  dynamic  stress  intensity  factors 
at  various  times,  obtained  by  the  dynamic  finite  element  analysis  [4],  are 
shown  in  Fig.  3.  As  seen  from  the  figure  the  profiles  are  nearly  linear  except 
very  near  the  crack-tip.  To  compare  with  the  static  case  which  does  not  include 
the  inertia  effect,  a  series  of  static  analyses  is  perforared,  with  the  boundary 
conditions  and  crack  lengths  corresponding  to  those  in  the  dynamic  analysis  [4] 
at  the  various  times.  Fig.  4  shows  the  crack  opening  profiles  in  the  static 
case.  Comparing  Figs.  3  and  4,  it  is  seen  that  the  crack  profiles  for  the 
bota  cases  .ire  very  similar,  but  with  different  amounts  of  crack  opening, 
at  the  respective  times. 

To  seek  i  correlation  between  the  dynamic  and  static  cases,  the  following 
eoor:icivnt  is  introduced  as 

-  -  (1) 

vii ere  :  is  the  crack-mouth  opening  (Fig.  1).  For  a  given  specimen  size,  this 
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coefficient  depends  on  not  only  the  crack  lc  "th  but  also  the  crack  speed  for 
the  dynamic  case*  while  the  coefficient  is  a  function  of  only  the  crack  length 
for  the  static  case.  The  variation  of  the  coefficients  for  both  the  static 
and  dynamic  cases  is  shown  in  Fig.  b.  An  excellent  linear  correlation  can  be 
seen  in  the  figure.  This  correlation  suggests  the  applicability  of  the  coef¬ 
ficient  obtained  by  the  static  analysis  to  the  evaluation  of  dynamic  stress 
intensity  factors. 

First  we  shall  seek  the  relation  between  the  static  stress  intensity  factor 
and  the  crack  mouth  opening  displacement.  According  to  Tada’s  handbook  [6], 
these  are  expressed  as 

K  •  */ira  F(a/W)  (2) 

L  2W 

<5  =  V(a/W)  (3) 

E’tT 

2 

where  E,x*E  for  plane  stress  and  E,=:E/(l-v  )  for  plane  strain.  From  the  above 
equations,  the  static  stress  intensity  factor  can  be  expressed  in  terms  of  the 
crack  mouth  opening  displacement  6: 

K  -  C  .  5  (4) 

I  s 

with 

cs  -  u(a/u>  <5> 

where 

bV  i/V)  *  (a/W)~**F(a/W)/V(a/W)  (6) 

The  function:*  F(a/W)  and  V(a/W)  have  been  reported  by  various  investigators. 
Among  them,  the  handbook  [b|  gives  the  tollowing  formulae: 
tor  S/'.J  -  4.0 

K(:x/V)  L. U‘)0  -  L. 735(a/W)  +  ».20(a/W)2  -  i4.13(a/W)3  +  14.57(a/W)4 

(0  a/W  <  0.6")  (i) 
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V(a/W)  =  0.76  -  2. 28(a/W)  +  3.87(a/W)2  -  2.04(a/W)3 


(l-a/W)‘ 


(any  a/W) 

and  for  pur  bending  (S/W=®); 


F(a/W) 


FT  a 

'Y  it  w 


_  ,ir  a. 
tan  (2  W> 


0.923+0. 199(l-sin(^  £)}4 
l  w 

,i ta. 

c°8(2W> 


(any  a/W)  (9) 

V(a/W)  =*  0.8  -  1.7 (a/W)  +  2.4(a/W)2  +  —  -y 

(1-a/W) 

(any  a/W)  (10) 

The  function  U(a/W)  can  be  calculated  by  substituting  Eqs.  (7)  and  (8), 
or  Eqs •  (9)  and  (10)  in  Eq.  (6).  Since  Eq.  (7)  is  invalid  for  a/W>0.6  and  the 
function  U(a/W)  is  less  sensitive  to  the  normalized  span  length  S/W,  the  formu¬ 
lae  for  pure  bending,  Eqs.  (9)  and  (10),  are  used  to  calculate  the  function 
U(a/W),  in  the  present  paper. 

The  variation  of  the  function  U(a/W)  with  the  normalized  crack  length  is 
shown  in  Fig.  6.  Dynamic  stress  intensity  factors  for  a  stationary  crack  under 
the  impact  loading  condition  can  be  estimated  from  Eq.  (4),  by  substituting  a 
value  of  U(a/W)  and  the  experimentally  measured  5. 

However,  for  a  dynamically  propagating  crack,  because  of  the  influence 
or  cue  crack  speed  on  the  displacement  field  near  the  crack-tip,  a  crack  speed 
correction  factor  is  required  to  determine  the  stress  intensity  factors.  The 
crack  opening  displacement  near  a  propagating  crack-tip  is  expressed  as  [71: 

V-)  -  ~tr  A(c)  (11) 

where  is  the  stress  intensity  factor  for  dynamically  propagating  crack, 
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G  is  the  shear  modulus,  r  is  the  distance  from  the  crack-tip,  and  A(c)  is  a 
dimens  ionless  function  of  the  crack  speed.  For  isotropic  materials,  the  factor 
A(c)  is  given  by 


A(c) 


a0162  -  u+e2)2 


»^l-(c/Cj)2;  j=l,2 


The  wave  velocities 


and  c^ 


are  given  by 


c 


1 


JS±1  6 

><-l  p 


and 


c 


2 


(12) 


(13) 


(14) 


(15) 


with 


3-4v 

(3-v)/(l+v) 


In  the  limit  when  c=0, 
A(0)  =  ^ 


for  plane  strain 
for  plane  stress 

the  value  of  the  function  is  given  by 


(16) 


(17) 


If  we  substitute  Eq.  (17)  into  Eq.  (11),  Eq.  (11)  becomes  identical  to  the 
crack  opening  displacement  for  a  stationary  crack  in  a  static  or  elastodyanamic 
field.  Comparing  both  the  stationary  and  propagating  cracks,  we  obtain  the 
following  relation: 


D 


A(0 1  .. 
Me)  NI 


(18) 


Thus,  :ao  i  tress  intensity  factor  for  a  propagating  crack  can  be  determined 
in  terms  <u  experimental ly  measured  *  as 


'ID 


i'  .  > 


;) 


(19) 


w  i  t  a 


)  - 


cn  «  U(a/W)  .  B(c) 

D  4^7 


(20) 


where  11(c)  is  the  crack,  speed  correction  l actor  and  given  by 

B(c)  =  A(0) ( A(c)  (21) 

The  variation  of  the  factor  B(c)  under  the  plane  strain  condition,  with  the 
normalized  crack  speed  (normalized  by  the  shear  wave  speed)  is  shown  in  Fig. 

7.  The  factor  B(c)  versus  c/c^  curves  depend  only  on  the  Poisson’s  ratio. 

The  value  of  B(c)  becomes  zero  at  c=c  where  c  is  the  Rayleigh  wave  speed. 

R  R 

If  the  crack  speed  is  relatively  slow,  i.e.,  c^O.lSc^,  a  value  of  1  can  be 
used  for  the  crack  speed  correction  factor  B(c)  allowing  an  error  within  2%. 

In  order  to  demonstrate  the  applicability  of  the  present  method,  the  crack 
mouth  opening  displacements  obtained  in  Ref.  [4]  are  used.  Fig.  7  shows  the 
variations  of  stress  intensity  factors  obtained  from  different  techniques. 

In  Fig.  7,  the  solid  line  shows  the  stress  intensity  factor  determined  directly 
as  a  variable  in  the  moving-singularity  element  procedure  [4],  and  the  dotted 
line  shows  that  calculated  by  Eq.  (19)  substituting  the  crack  mouth  opening 
displacements  determined  in  Ref.  [4].  As  demonstrated  in  Ref.  [3],  the  moving- 
singularity  finite  element  method  gives  very  accurate  dynamic  stress  intensity 
factors.  As  seen  from  the  figure  the  present  estimation  procedure,  Eqs.  (19,  20), 
yields  results  which  agree  excellently  witli  those  of  the  moving-singularity 
element  procedurc[4]  when  the  crack  is  stationary.  When  the  crack  propagates, 
however,  the  overall  variation  of  the  presently  estimated  result  is  good,  although 
the  estimated  result  appears  to  oscillate  slightiv  around  the  result  of  the 
rnovin*;  singularity  element  procedure  [4).  For  a  given  material,  it  is  noted 
that  the  present  s imple-es t imat ion  method  gives  better  results  for  a  shorter  crack. 

or  otuer  types  of  specimen  which  have  a  relatively  short  edge  crack  with- 
Mt  i nv  loading  on  the  crack  surface,  the  present  method  can  also  be  applied 


to  determine  the  dynamic  stress  intensity  factors  from  the  mouth  opening 
displacements  of  the  propagating  crack. 
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Figure  Captions: 

Fig.  1:  Dynamic  tear  test  specimen 

Fig.  2:  Crack  growth  history 

Fig.  3:  Crack  opening  profiles  in  the  dynamic  case 

Fig.  4:  Crack  opening  profiles  in  the  static  case 

Fig.  5:  Variation  of  the  coefficients  C(=Kj./6) 

in  the  dynamic  and  static  cases 

Fig.  6:  Function  of  crack  length  U(a/W) 

Fig.  7:  Crack  speed  correction  factor  B(c) 

Fig.  8:  Comparison  of  dynamic  stress  intensity  factors 
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Thickness  B  =  l5.8mm  Mass  Density  p  -  7840  Kg/m 


(  )  =  Ki  (MNm'1'5) 


(32.6)  a/w  =0.25  (25  /isec) 


(65.2)  a/w=0.25  (50  Msec) 


(97.8)  a/w=0.25  (75  Msec) 


(1334)  a/w= 0.299  (100  Msec) 


a/w  =  0.546  (125/isec) 
(153.4) 


a /w=  0.763 

(150  asec) 
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